Draft version October 26, 2012 

Preprint typeset using L£-T^X style cmulatcapj v. 12/16/11 



A MULTI-SPACECRAFT VIEW OF A GIANT FILAMENT ERUPTION DURING 26/27 SEPTEMBER 2009 

Sanjay Gosain 1 , Brigitte Schmieder 2 , Guy Artzner 3 , Sergei Bogachev 4 , and Tibor Torok 6 

Draft version October 26, 2012 

ABSTRACT 

We analyze multi-spacecraft observations of a giant filament eruption that occurred during 26 and 
27 September 2009. The filament eruption was associated with a relatively slow coronal mass ejection 
(CME). The filament consisted of a large and a small part, both parts erupted nearly simultane- 
ously. Here we focus on the eruption associated with the larger part of the filament. The STEREO 
satellites were separated by about 117° during this event, so wc additionally used SoHO/EIT and 
CORONAS/TESIS observations as a third eye (Earth view) to aid our measurements. We measure 
the plane-of-sky trajectory of the filament as seen from STEREO-A and TESIS view-points. Us- 
ing a simple trigonometric relation, we then use these measurements to estimate the true direction 
of propagation of the filament which allows us to derive the true i?/i?Q-time profile of the filament 
apex. Furthermore, wc develop a new tomographic method that can potentially provide a more robust 
three-dimensional reconstruction by exploiting multiple simultaneous views. We apply this method 
also to investigate the 3D evolution of the top part of filament. We expect this method to be useful 
when SDO and STEREO observations are combined. We then analyze the kinematics of the eruptive 
filament during its rapid acceleration phase by fitting different functional forms to the height-time 
data derived from the two methods. Wc find that, for both methods, an exponential function fits the 
rise profile of the filament slightly better than parabolic or cubic functions. Finally, we confront these 
results with the predictions of theoretical eruption models. 
Subject headings: CME, Filament eruption 



1. INTRODUCTION 

Coronal mass ejections (CMEs) are huge expulsions 
of plasma and magnetic field from the solar corona into 
interplanetary space. They are often accompanied by 
the eruption of a filament or prominence, which becomes 
visible as the core of the CME in coronagraph observa- 
tions, and by a flare that occurs almost simultaneously 
with the eruption. It is now well accepted that these 
three phenomena are different observational manifesta- 
tions of a more general process, nam ely a local d isrup- 
tion of the coronal mag netic field (e.g.. lForbesl2000f ). The 
detailed mechanisms that trigger and drive such disrup- 
tions are, however, still controversial, and a large num- 
ber of theoretical models has been put forward in the 
past decades (for recent reviews, see lAmari fc KbMM . 
lAulanier et al1[20i0l : !Forbesll2lnoh . 

Early observations indicated that there are two distinct 
classes of CMEs, namely fast (or impulsive) ones, origi- 
nating in active regions and associated with flares, and 
slow (or gradual) ones, associated with large prominence 
eruptions outsi de of active regions and no, or no sig- 
nifica nt, flaring {MacQuecn fc Fisherlll983HSheelev et al.l 
Il999h . Consequently, it has been suggested that differ- 
ent eruption mechanisms may be at work in these two 
types of eruptions. However, the analysis of considerably 
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larger data sets in the SoHO era revealed a continuous 
distribution of CME velocities with a single peak (e.g., 
I Zhang fc Derdl2~006[ ). indicating that both fast and slow 
CMEs are driven by the same physical mcchanism(s). 
This is supported by the considerable range of CME kine- 
matics that could be modeled based on a single physica l 
mechanism (jChen fc Kralll 120031 ITorok fc Kliernl l2007h , 
as well as by the fact that large prominence eruptions 
outside active regions can produce loops and ribbons that 
are morphologically similar to those seen in flare-related 
CMEs. The majority of prominence-related CMEs are 
most likely not associated with flares simply because the 
magnetic fields in their source regions are too weak to 
produce significa nt emission i n Ha and in EUV wave- 
lengths (see, e.g.. lForbesll200Q[) . 

Virtually all theoretical models describe CMEs as coro- 
nal magnetic flux rop es that are anchore d in the dense 
photosphere (see, e.g.. lGibson et al.ll2006f ). although it is 
debated whether a flux rope is present in the corona prior 
to an eruption or is formed during the eruption process. 
The expulsion of a flux rope into interplanetary space as a 
CME has been explained by, e.g., the continuous increase 
of poloidal flux in the rope due to flux injection from 
the c onvection zone into the corona (e.g.. iChen fc KraH 
I2003h . by i deal MHD instabilities like the helical kink 
instability ()Fanll2005l; ITorok fc Klieml l2005D or the torus 
instability ( Kliem fc To rok 2006), or by the combination 
of a "loss of equilibrium" of a flux rop e and magnetic re- 
conn ection occurring in its wake (e.g. JForbes fc Isenbergl 
Il991h . Other models invoke reconnection from the begin- 
ning of the eruption, as for example the "tether cutting" 
(e.g..lMoore et aLll2001f) and "magnetic breakout" (e.g., 
lAntiochos et al.lll999f ) models. 

Using observations to support or reject specific models 
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for a particular eruption is difficult for several reasons. 
First, many models employ more than one physical mech- 
anism, resulting; in a partial overlap between them (see 
lAulanier et al.1 l2010f ) . Second, several distinct mecha- 
nisms may occur almost simultaneously in an eruption, 
in particular in complex events, making it difficult to 
estab lish which one is the main driver of the eruption 
(e.g.. iWilliams et aI]|2005D . Third, the models predict a 
very similar evolution for the main phase of an eruption, 
i.e., for the evolution after the impulsive flare phase and 
initial rapid acceleration of the ejecta. When looking 
for clues for possible eruption mechanisms, one therefore 
often focuses on the early eruption phase, for example 
on the morpholog y and timing of pre- fl are Ha and EUV 
brightening (e.g.. [Chandra et al.ll2010T) . For eruptions 
associated with large quiescent prominences, as the one 
studied in this paper, such signatures are, however, often 
not available. 

Another possibility to obtain information about the 
mechanisms at work in an eruption is to study its 
kinematic properties, in particular the early rise phase 
(jSchriiver et al.l 120081 ). Eruptions typically start with 
the slow rise of a filament or prominence and/or over- 
lying loops, at an approximately constant velocity of 
a few km/s, which is followed by the rapid accelera- 
tion of the ejecta to several 100 km/s. The accelera- 
tion should initially follow some functional dependence, 
but will t hen saturate and decr ease afterwards (see, e.g., 
Fig. 3 in IGallagher et al.l 12001 . Profiles of the initial 
acceleration phase, if extractable from measured height- 
time data with sufficient coverage, can be compared with 
predictions of theoretical eruption models. 

Observed rise profiles of the early phase of fil- 
ament eruptions and CMEs ha ve been fitted by 
constant-acceleration curves (e.g.. IGilbert et al.l |2000; 



iKundu et al1l200l iChifor et al.ll2006D. power-law. h (t) tx 
t m , with 3.0 5, m < 3.7 (e.gTlAlexander et all [20021: 
iSchriiver et all 120081 [Liu et al.1 20091) . and exponentia l 
functions ([Gallagher et al.l 120031: IWilliams et ail I2005D . 
As for theoretical models, the functional dependence of 
eruption trajectories has not often been reported yet, and 
systematic investigations of its parameter dependence 
are quite rare so far. Still, some model predictions can be 
inferred from the literature. For example, an exponential 
early rise is naturally expected if an erupti on is initially 
drive n by an ideal MHD instability (e.g., ITorok et al.l 
l200l . A power-law dependence with m = 2.5 has been 
found for the trajectory in a two-dime nsional (2D) ver- 
sion of the loss of equilibrium model ([Priest fc Forbes! 
12002ft . and a parabolic rise (i.e., constant acceleration) 
was reported for a simulation o f the breakout mode l 
([Lynch et al.ll2004f ). Furthermore. ISchriiver et al.1 ([20081 ) 
showed that a velocity perturbation at the onset of the 
rapid acceleration of an eruption can change the resulting 
functional dependence of the trajectory (from exponen- 
tial to near-cubic for the case of the torus instability they 
studied), given that this perturbation is sufficiently large 
(somewhat larger than typical velocities observed during 
the initial slow rise p hase of an eruption) . We refer to 
ISchriiver etUI ([20081 ) for further details. 

Prominence eruptions and CMEs have been observed 
for a long time with various ground and space-based 
instrument. For exam ple, the LASCO coronagraphs 
(|Brueckner et al.l [1995) onboard the SoHO spacecraft 



have been observing thousands of CMEs in white-light. 
However, one of the limitations of LASCO, and of other 
instruments that are located at, or close to, Earth, is that 
they can only obtain 2D observations, projected onto the 
plane-of-sky (POS). Height-time data of eruptions are 
particularly hampered by this, since obtaining the cor- 
rect radial rise velocities requires the knowledge of the 
true 3D trajectories. 

STEREO mission ([Kaiser et al.l 120081) was launched 
during solar minimum and theref ore initial studies 
were of quiescent fil a ment eruptions dGosain et al.l[2009t 
lArtzner et all MHI IGosain fc Schmieden I2010TP With 
STEREO 3D reconstructions, it is, in principle, possi- 
ble to derive the shape of the prominence, its twist or 
writhe, and its true trajectory when it erupts. These 
properties can b e useful for comp a risons with model 
predictions (e.g ITorok et al.l 120101 : iKliem et al.l 120121 : 
iZuccarello et al.ll2012f ).~ Other case studies of quiescent 
fil ament erup t ions o bserv ed with STEREO a re reviewed 
in iBemporadl (1201 If) and lAschwandenl ([20111 ) and those 
of CMEs in lMierla et al.l (|2010l) . 

The main difficulty in stereoscopic reconstruction 
arises when the separation angle between the two 
STEREO satellites is large, because then it becomes dif- 
ficult to iden tify the same feature i n both views un- 
ambiguously ([Thompson et all l2012t) . In such cases 
one has to use complementary ob servations from other 
instru ments, such as SoH O/EIT dDelaboudiniere et al.1 
fl995l) . CO RONAS/TESIS (iKuTin et al.l 120091) and now 
SDO/AIA (jLemen et all l2012t ). The addition of the 
Earth view to STEREO views makes life easier as it pro- 
vides us (a) more than one stereoscopic pair, and (b) 
smaller separation angles. Also, in a special circumstance 
when a structure (filament/promincence) is not visible in 
one of the STEREO satellites (for e.g., remains hidden 
behind the limb during initial rising phases), adding an 
Earth view to one of the STEREO satellites allows us to 
make stereoscopic reconstruction. So me recent exam ples 
where three views have been used are lLi et all (1201 ID and 
iFeng et al.l (pOlU ). 

In this paper we will present such an example. We 
present He II 304 A observations of a large filament erup- 
tion that occurred during 26-27 September 2009. The 
observations were taken from the twin STEREO satel- 
lites and were combined with complementary observa- 
tions from SoHO/EIT and CORONAS/TESIS, giving us 
a third view i.e., the Earth view of the event. The fila- 
ment eruption was seen as a limb-event by STEREO-A, 
EIT, and TESIS, while it was seen as an on-disk event 
by STEREO-B. 

Based on these multi-spacecraft data, we derive the 
R/Rq or simply height-time profile of the erupting 
prominence by different methods. First, we derive in- 
dependently the POS height-time profiles of the promi- 
nence top as viewed from STEREO-A and TESIS. Sec- 
ond, we apply a trigonometric relation to simultane- 
ous STEREO-A and EIT observation of the prominence 
and estimate the propagation direction of the filament. 
Knowing the propagation direction we can derive true 
height-time profile of the filament top. Finally, we apply 
a three-dimensional stereoscopic reconstruction method 
based on Marinus projections to derive the height-time 
profile, we call this method "tomographic method" as it 
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can use simultaneous multiple views for reconstruction. 
We then fit these height-time profiles by functional forms, 
viz., parabolic, exponential and cubic and compare the 
results. 

The paper is organized as follows: Section 2 sum- 
marizes the observations. In Sections 3 and 4 we de- 
scribe the trigonometric and tomographic methods, re- 
spectively results derived from these methods. In section 
5 we discuss the results as well as the potential of the 
newly developed tomographic method. 

2. OBSERVATIONS AND METHOD OF DATA ANALYSIS 
2.1. He II 304 h Filament Observations 

During 26-27 September 2009 a large filament eruption 
was observed near the north-east solar limb. The obser- 
vations were obtained in the He II 304 A wavelength at 
a cadence of ten minutes by the SECCHI/EUVI instru- 
ment aboard the STEREO-A and B satellites. Figure Q] 
(top panel) shows on-disk observations by STEREO-B 
during the early phase of the eruption. Two filaments 
can be seen, a long one (LF) located at 23-34° N, and a 
short one (SF) located at 18-32° N. The bottom panel 
shows a reconstruction of the photosphcric magnetic field 
distribution as it would have been seen from STEREO-B 
(at 19:17 UT on 26 September 2009, using B a = 5.54° 
and CML = 164.4°, where B and CML are the Car- 
rington coordinates of the disk center in the STEREO-B 
view) . It can be seen that both filaments follow the same 
polarity inversion line (PIL), suggesting that they were 
both located within an extended filament channel. The 
magnetic field on cither side of the PIL was weak (+/-3 
Gauss). In the limb view of STEREO-A, LF starts to rise 
earlier than SF (see movie 1), while both eruptions seem 
to occur simultaneously in the disk view of STEREO-B 
(movie 2). Since LF is much more prominent than SF, 
we focus in thi s paper on the evolution of LF and refer to 
Li ct al. (2010) for further details on the eruption of SF. 
In the following, "filament" or "prominence" therefore 
refers only to LF. 

During this event the STEREO satellites were sepa- 
rated by about 117° and the main part of the filament 
was hidden behind the solar limb in the STEREO-A 
viewpoint until about 18:20 UT on 26 September 2009. 
A 3-D stereoscopic reconstruction in combination with 
STEREO-A before 18:20 UT was therefore not possible 
from either STEREO-B or Earth view (EIT and TESIS). 
He II 304 A images from SoHO/EIT at a cadence of 
six hours and CORONAS/TESIS images at a cadence 
of ten minutes provide the Earth view, which are used 
for 3D reconstruction with STEREO-A after 18:20 UT 
Sept 26 using the tomographic method described in sec- 
tion 3.3. The He II 304 A images of the filament from 
the four instruments are shown in Figure [2] The Earth- 
Sun-STEREO-A and Earth-Sun-STEREO-B angles on 
26 September 2009 were 61° and 56°, respectively. The 
observations of the filament top from two different van- 
tage points (TESIS and STEREO-A) allows us to observe 
the POS evolution in a piecewise continuous manner from 
TESIS (00:00 to 22:00 UT) and STEREO-A (19:10 to 
23:10 UT), as shown in Figure M (top left panel). A 
simple trigonometric method, described in section 3.2, 
is used to triangulate the true propagation direction of 
the filament apex. Knowing this angle the POS height- 



time profiles are corrected to derive the true height-time 
profiles as shown in the Figure [6] (bottom left panel). 

The filament LF (as seen from STEREO-A; Figure [3]) 
suggests a shee t-like mo r pholo gy. The stereoscopic re- 
construction bv lLi et all ()2010l ) (see their Figure 2) also 
infers a sheet like structure. We outline its appar- 
ent edges by dashed (red) and dotted (yellow) lines in 
STEREO-A and B views. A careful inspection of the 
legs of the prominence in STEREO-A images suggests a 
twisted morphology of the legs. However, a quantifica- 
tion of the twist of the filament sheet is not possible in 
the present case. 

2.2. CME observations 

STEREO-A observed the CME associated with the fil- 
ament eruption on 26/27 September 2009 with its coron- 
agraphs, COR1 (1.4 to 4 R ) and COR2 (2.5 to 15 R ). 
The CME was not seen by STEREO-B since it was di- 
rected towards it and was perhaps too faint to be seen 
as a halo CME. The LASCO/C2 (1.5 to 6 R Q ) and C3 
(3.7 to 30 Rq) coronagraphs also observed the CME. The 
time of arrival of CME in C2 was 23:06 UT on Sept 26, 
in C3 was 14:18 UT on Sept 27, and in HI1 was on 21:29 
UT on Sept 27. 

The propagation angle of the CME leading edge is de- 
rived using same procedure as described in section 3.2 
and applied to filament apex. The projection correc- 
tion to the POS height-time profiles of the CME (top 
right panel of Figure EJ is then applied to derive the true 
height-time profiles (bottom right panel of Figure |6]) . 

3. HEIGHT-TIME PROFILE OF ERUPTING FILAMENT 
AND CME 

3.1. Plane- of- sky Measurements 

We make plane-of-sky (POS) measurements of the fil- 
ament apex (marked by '+' symbol in the panels of Fig- 
ure [2]) using TESIS and STEREO-A observations. The 
Rpos / Rq profile of the filament apex measured from 
the two views (from STEREO-A and TESIS) is plotted 
in the top left panel of the Figure |6] Similarly, the top 
right panel of the Figure [5] shows the Rpos/Rq profile of 
the leading edge of the CME measured with STEREO-A 
and LASCO coronagraph observations. The difference 
in the Rpos/Rq profiles of the filament apex and CME 
leading edge as seen by different satellites is apparent, 
since the measurements correspond to two different van- 
tage points. 

3.2. Estimation of the True Height-Time Profile using 
Simple Triangulation 

The POS measurements, Rpos/Rq, of the filament 
apex and the CME leading edge, described above, can 
be corrected for projection effect if we know the angle 
between the real trajectory of the erupting feature and 
the POS in the observer's frame of reference. The true 
R/Rq profile is related to the Rpos/Rq profile, mea- 
sured in the POS, by (R/R Q )cos 9 = Rpos/Rq, where 
9 is the angle between the real trajectory and the POS, 
referred to as propagation angle henceforth. 

Here we apply a simple trigonometric relation, using 
image pairs from STEREO-A and TESIS/EIT (hence- 
forth, Earth View or EV), to estimate the propagation 
angle. We then use this information to derive the true 
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R/Rq profile of the erupting filament and the CME. A 
simple assumption made here is that the propagation an- 
gle remains unchanged during the time of the measure- 
ments. We verified this assumption by computing the 
propagation angle using STEREO-A and TESIS pair at 
later times and found that the angle remains the same 
(see Table 1). 

We explain the trigonometric procedure here briefly. 
The illustration in Figure [4] shows the geometric set- 
ting of the two STEREO satellites and the EV with 
respect to the filament. The projected height of the 
top part of the filament is h a and hb in the POS of 
STEREO-A and EV, respectively. S is the separation 
angle between STEREO-A and EV (61 degrees), and 
a and (3 are the angles that the top of the filament 
apex (trajectory) makes with respect to the POS. We 
obtain the angles using the relations S = a + (3 and 
h a /hb = cos(a)/cos(S f — a). Knowing these angles, we 
can apply corrections to the POS heights h a and hb to 
obtain the true height hxrue = h a /cos(a) = ft.(,/cos(/3). 

Figure [5] shows an example of two stereoscopic image 
pairs, i.e., STEREO-A (left panel) and EIT (right panel), 
observed almost simultaneously. The two images are in 
cpipolar view. The segments h a and hb measure the top 
part of the filament as viewed from two vantage points. 
Knowing h a , hb, and S, we determine a and /3 to be 41 
and 20 degrees, respectively. Since only one stereoscopic 
pair is available between SOHO/EIT and STEREO-A, 
we make use of TESIS data to make pairs with STEREO- 
A. The observations of STEREO-A and TESIS are not 
synchronized in time and both instruments follow a dif- 
ferent time cadence (see Table 1). Assuming that the fil- 
ament did not evolve significantly within the small time 
differences, we could make 8 near-simultaneous stereo- 
scopic pairs of TESIS and STEREO-A. The timings of 
these pairs and the value of angles a and /3 deduced using 
these pairs is given in Table 1. It may be noticed that 
during the observed time interval the propagation angle 
does not change and therefore it is possible to correct 
the observed POS height-time profiles for propagation 
angle using single value of a and /?. In section 3.3.3 we 
show that these values are consistent with other 3D re- 
construction methods. The corrected height-time profile 
of the filament apex is shown in the lower left panel of 
the Figure [6] 

Similarly, applying this method to the CME leading 
edge we deduce angles a and (3 to be 36 and 25 degrees, 
respectively. The true height-time profile of CME lead- 
ing edge after correcting for these angles is shown in the 
lower right panel of the Figure [5] It is interesting to note 
that the direction of propagation of the filament apex and 
the CME leading edge differs by about 5 degrees. The 
filament typically forms a core in the three part CME 
structure. However, since the CME leading edge is more 
extended, i.e., the front surface of a tear-drop shaped 
bubble in which filament forms a trailing part, the differ- 
ence of 5 degrees is small considering the large angular 
extent of the filament and the associated CME. Another, 
interesting point about this method is that the two curves 
merge into one (as seen in the combined curves in lower 
panels) only for a unique pair of a and f3 angles, where 
the sum of the two angles (S — a + /3, in this case equal 
to separation between Earth and STEREO-A, i.e., 61.5 
degrees) is well constrained by the known separation an- 



gle between the two vantage points. For any other pair 
of these angles the two curves did not merge into one. 
Thus, just by knowing the separation angle between two 
vantage points and the respective POS height-time pro- 
files, one can iteratively adjust the angles (in fact only 
one of the angles, as the two angles a and (3 arc sim- 
ply, a and S — a), until the two curves merge as one. 
This procedure also gives the same solution for a and (3. 
These true height-time profiles are then used for deriving 
the velocity and acceleration profiles of the filament and 
CME, which is described in the following sections. 

3.2.1. Estimating the duration of Rapid acceleration phase: 

In this section, we use the true height-time curves, 
shown in the lower panels of Figure [6] to derive the ve- 
locity, acceleration and jerk (rate of change of acceler- 
ation, following Schrijver et al. (2008)) profile of the 
filament and CME. The latter can then be compared to 
the predictions of theoretical eruption models described 
in the Introduction. Since the acceleration and jerk arc 
the higher order time derivative of the trajectory, er- 
rors in the measured data amplify strongly, so one is 
typically forced to smooth the data before calculating 
acceleration curves, for example using spline smoothing 
(jVrsnak et all [20071) . Here we will use a different ap- 
proach: we first fit a fourth order polynomial of the form 
H(t) = a + bt + ct 2 + dt 3 + et 4 to the height-time data. 
We then use this smooth curve to obtain the velocity, 
acceleration and jerk profiles. These profiles are shown 
in Figure [7| for the filament and the CME leading edge 
in the left and right columns, respectively. 

It is to be noted that the underlying physical mech- 
anism responsible for the eruption determines the func- 
tional form of the height-time profile only initially, i.e., 
during the phase when acceleration is growing, but not 
yet saturating. Once the acceleration starts to saturate, 
the functional form is changing. From the FigurcJT] we 
notice that the acceleration profiles as estimated from 
the 4th order polynomial are quite different for the fila- 
ment and the CME. While the acceleration of the CME 
leading edge is higher than that of the filament, the rate 
of change of acceleration, i.e, the value of jerk for the 
CME is declining, in contrast to the filament. This sug- 
gests that during the time interval of the CME data the 
increase of the acceleration of the leading edge is slowing 
down, while the growth of the acceleration of the fila- 
ment is still increasing. We therefore restrict our fits of 
different functional forms to the filament only and not to 
the CME. 

3.2.2. Fitting functional forms to the Filament rapid 
acceleration phase: 

Before we fit the functional forms we make an estima- 
tion of the optimum time interval which corresponds to 
rapid acceleration phase of the filament. To get the first 
estimation we took the start time where the acceleration 
starts to grow from zero and the end time as the last 
data point. We then fine-tune our estimation of the time 
interval of the rapid acceleration phase by varying their 
start and end times and observing the resulting quality 
of the overall fits of all three functional forms. The inter- 
val leading to an overall best fit quality for all functional 
forms is marked by vertical dashed lines in Figure [71 
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Within this rapid acceleration interval we then find 
the best fit functional form as described below. We 
fit the three different functional forms: (i) parabolic, 
H(t) = a + bt + ct 2 , (ii) exponential, H(t) = ae bt +c, and 
(iii) cubic, H(t) = a + bt + ct 3 to the height-time profile 
during the rapid acceleration phase. The left panel of 
the Figure M shows these fits. The reduced chi-square \ v 
values arc shown at the top left corner of each panel. The 
weights that we apply to the data points for fitting arc 
taken to be W = 1/cr 2 , where a is the standard deviation 
of the measurement error. We assume a Gaussian distri- 
bution for the latter. The pixel size of STEREO/EUVI 
is about 2 Mm. For the filament (or prominence), which 
is typically quite diffuse in He II 304 A images, we con- 
sider 3 pixels, i.e. 6 Mm as the 1-tr error. The actual 
errors may be somewhat different. However, while us- 
ing different values of 1-cr will lead to different values of 
Xv, the relative values of \ v between different functional 
forms will remain the same. The fits of the three func- 
tional forms shown in the left panels of Figure |H1 clearly 
favor an exponential rise of the filament, with a relatively 
better value of Xv 

3.3. The 3-D reconstruction by using Marinus 
projection 

3.3.1. The Method 

Here we describe a new tomographic method for the 
3D reconstruction. We used simultaneous views of the 
filament from STEREO-A and B and TESIS in He II 304 
A wavelength. The essence of the method is as follows. 

A continuum intensity image of the sun, I(x,y), can 
be easily projected into heliographic coordinates 1(1, b). 
This projection is also known as equidistant cylindrical 
or Marinus projection. Since the continuum intensity 
I(x, y) corresponds to the solar photosphere, each point 
(x, y) on the intensity image can be associated with heli- 
ographic coordinates (/, 6), assuming a spherical sun with 
radius, R = Rq. A common feature on the solar disk, 
such as a sunspot should then correspond to the same 
Carrington latitude-longitude, no matter what the view- 
ing angle of the sun is. However, for coronal images like 
in He II 304 A the intensity features corresponding to 
filaments, spicules etc. do not lie on the same sphere but 
are elevated structures in 3D. Thus, a common feature 
like a filament or coronal loop will correspond to location 
(7i, b\) and (I2, ^2) in heliographic projection of the coro- 
nal images obtained from different viewing points 1 and 
2, respectively (when R = Rq is assumed). Conversely, 
if the heliographic projection is attempted assuming the 
Sun to be a sphere of radius larger than one solar radius 
and a correct radius of the sphere is assumed (equal to 
the altitude of the feature) then the we should get l\ = h 
and b\ = bi for the common feature. 

Thus, generating the generalized Carrington maps for 
different assumed radii of the spherical grid, using a 5 
Mm step from i?=700 Mm to i?=1500 Mm, and compar- 
ing the Carrington coordinates (latitude-longitude) of a 
recognizable common feature, such as filament apex, in 
the three Carrington maps (one for each viewing angle) 
until they all agree gives us a solution for the 3-D co- 
ordinates of the feature. We found that this step size 
of 5 Mm gives an optimum choice to arrive at the best 
agreement for the generalized coordinates of a recogniz- 



able feature. Thus the generic accuracy of the method 
can be assumed to be about 5 Mm. 

We geometrically consider both intersections of the line 
of sight with the reference sphere. When the radius of 
the reference sphere is equal to the chromospheric radius, 
we take into account the single point located physically 
in front of the POS. When the radius of the reference 
sphere is greater than the chromospheric radius, we have 
in principle to take into account both intersection points 
located respectively in front of and behind the POS. That 
is why in the top panel of Figure 8 the prominence, pro- 
jected on the far side of the reference sphere, behind the 
POS of STEREO-A, appears as reversed from right to 
left with respect to the direct view in Figure 2. In addi- 
tion to that, the far side of the solar disk appears as a 
dark, missing disk in top panel of the Figure |SJ 

3.3.2. Advantages and Limitations of the Method 

It is well known that all stereoscopic reconstruction 
methods are limited by the ambiguity in recognizing a 
common feature in different views. Further, a common 
limitation that arises with any 3D stereoscopic recon- 
struction technique is when the apex point from two 
viewing angles may be different. Such situations would 
lead to a systematic error in the reconstructed 3D coordi- 
nates. However, we expect such errors to be less severe in 
our case because (a) the filament studied here has a large 
extension in longitude which is rising globally as a whole, 
so height-time profile of several neighboring points along 
the filament will be similar, and (b) using combination of 
STEREO-A and TESIS (separation angle 61°) as com- 
pared to STEREO-A and STEREO-B (separation angle 
117°), we reduce the errors. Although, such systematic 
errors cannot be avoided, the time derivative of measured 
altitude and hence the derived velocity and acceleration 
should not be affected severely as long as the system- 
atic error remains similar in magnitude. Therefore, for 
studying the kinematic evolution of erupting prominence 
such reconstruction methods may still be applicable, with 
aforesaid limitations. 

The Figure [SJ shows selected parts of the three gen- 
eralized Carrington maps corresponding to STEREO- 
A, TESIS and STEREO-B views, generated assuming 
radii of the spherical grid to be 1245 Mm, i.e., 545 Mm 
above the solar surface. At the choice of this radius 
the common feature, i.e., the filament apex marked by a 
square box, corresponds to the same Carrington latitude- 
longitude coordinates in the different views. In principle, 
two stereoscopic views are sufficient for the application 
of this method. However, adding more views increases 
redundancy (for example in present case more emphasis 
is given to TESIS and STEREO-A for constraining re- 
construction) and therefore may add to its robustness. 
In future, we plan to apply this method to the events 
observed simultaneously by the two STEREO /EUVI in- 
struments along with the high- resolution SDO/AIA ob- 
servations. 

3.3.3. Comparison with SCC_MEASURE and Simple 
Triangulation Method 

The 3D reconstruction of the fi lament s t udied in this 
paper, was also carried out by iLi et al.l (|2010ft using 
SCC_MEASURE procedure (developed by W. Thomp- 
son). They used STEREO-A and STEREO-B pair for 
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their reconstruction. Further, they reconstructed many 
(12) points along the filament body (their Figure 5). 
Their points 6,7 and 8 correspond to the top of the fil- 
ament and one can notice that the height evolution of 
these points (their Figure 6(a)) is quite similar to each 
other (within ±6Mm), though the points are separated 
spatially, this is due to the large scale uniform evolution 
of the filament. For comparison, we overlay the altitude 
data p oints of location 7, as reconstructed by iLi et al.l 
(|2010f) in our height-time plot shown with red colored 
symbols in the top-right panel of Figure 9. The re- 
construction from two independent methods agrees quite 
well, considering the general scatter in the reconstructed 
coordinates. 

On the other hand, a poorer match is expected be- 
tween the true 3D reconstruction methods and the sim- 
ple triangulation method since the latter only estimates 
the propagation angle and not the 3D coordinates of the 
filament. The height derived from simple triangulation 
method shows a systematic offset with respect to the true 
height derived from 3D reconstruction methods. Apart 
from the systematic offset, the profile of the derived speed 
and acceleration should however remain unaffected, as 
these depend upon the shape of the curve. This is evi- 
denced in a similar fit quality of the height-time profile 
by both the methods to different functional forms (see 
section 3.3.4). Also, the method is straightforward and 
relies in tracking a common feature in the images taken 
from same vantage point albeit at different times. The 
natural advantage is that it is easy to track a common 
feature in time if the time difference between two images 
is not very large. 

Further, it may be noted from Table 1 that the angles 
a and /3 are not changing s ignificantly. In the 3D recon- 
struction by (|Li et al.ll201Q[ ) (their figure 4, right panels) 
it can be seen that a propagation angle of ^20° in front 
of the East solar limb is deduced and is not changing 
significantly. Also, the Carrington longitude of the fila- 
ment apex reconstructed using Marinus method (Table 
1.) shows a small variation in longitude of ~ 4°, while 
the mean value of the longitude, ~ 150°, corresponds 
to an angle of ~ 20° in front of the East solar limb, in 
agreement with lLi et al.l (|2010f ) and angle /3 from simple 
triangulation method. 

3.3.4. Rapid acceleration phase and its functional form: 

We use this tomographic 3-D reconstruction method 
based on Marinus projection to obtain the 3-D trajec- 
tory of the filament apex. The rise of the altitude of the 
filament apex is fitted for different functional forms. The 
time interval of rapid acceleration phase is taken to be the 
same as estimated in section 3.2.1. The height-time curve 
and the fitted parabolic, exponential and cubic functions 
to it are shown in the right panel of Figure [9] from top to 
bottom, respectively. The reduced chi-square \ v values 
are shown at the top left corner of each panel. It is found 
that an exponential form fits the observations relatively 
better as compared to the other functions. The exponen- 
tial function was also found to fit the rapid acceleration 
phase relatively better than other functions in section 
3.2.2, where simple triangulation method was used. This 
is shown in the left panel of Figure [9] 

4. DISCUSSION AND CONCLUSIONS 



In this paper we analyzed the observations of a large 
erupting quiescent filament which was observed from 
three vantage points by STEREO-A, B and the EV 
(SoHO/EIT and TESIS). The filament rose slowly for 
several hours before accelerating rapidly and erupting in 
two parts, a large and a small filament. We analyzed 
the kinematics of the large filament, whose true trajec- 
tory was derived by two methods: one simple triangula- 
tion method and another newly introduced tomography 
method. The new tomographic method can potentially 
take advantage of simultaneous observations from multi- 
ple vantage points to constrain the reconstructions bet- 
ter. After deriving the true trajectory by the two meth- 
ods we fitted the height-time curves with different func- 
tional forms and compared the results with predictions 
of theoretical eruption models. 

The key points in the observational analysis can be 
summarized as follows: 

1. The eruption involved two filaments, a large one 
and a small one, which were located above the same 
polarity inversion line, suggesting that they were 
embedded in the same, elongated filament channel. 
The photosphcric magnetic field strengths at the 
location of the filaments were weak (up to about 3 
Gauss). The two filaments erupted almost simulta- 
neously. In the present analysis we focused on the 
eruption of the more prominent large filament. 

2. We used two different approaches to derive the true 
height-time profile of the filament. First we used a 
simple triangulation method to determine the an- 
gle which the filament trajectory makes with re- 
spect to the planc-of-sky (POS) and applied cor- 
rection to the plane-of-sky height-time profile to 
derive true height-time profile. Second, we used to- 
mographic approach where we make Marinus pro- 
jections of the three views of the sun on spheres of 
radii larger than the solar radii so as to arrive at 
a common latitude-longitude position of a common 
feature (filament apex) in all maps. The advantage 
of the first method is that once we know the prop- 
agation angle with respect to POS from triangula- 
tion, we can go back and forth in time and correct 
the POS height-time profile obtained with even one 
satellite, i.e., durations when only one view is avail- 
able, for e.g. when in one of the stereoscopic pairs 
the filament is behind the limb or out of the FOV. 
However, the method assumes that the propagation 
angle of the filament with respect to POS does not 
change substantially over the time of observations. 

3. During its early rise phase, the filament exhibits 
the morphology of a twisted sheet. However, its 
chirality could not be inferred from the images. 

4. We derived the acceleration and jerk (rate of 
change of acceleration) profiles for the filament and 
the CME (Figure EJ. It is believed that the ini- 
tial rapid acceleration phase, when acceleration is 
growing, may be suggestive of the physical mech- 
anism behind the eruption (Schrijvcr et al. 2008). 
However, the acceleration curve must be growing 
and not saturating or slowing down, in other words 
the jerk should be increasing. By studying the jerk 
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profiles in Figure [7] for the filament and the CME 
we decided to fit different functional forms to the 
filament observations only and not to the CME, be- 
cause the jerk profile of the CME suggests that its 
acceleration is already saturating. Since the CME 
observations are available only when it enters the 
coronagraph's FOV, which is much later than the 
observations of the filament eruption, we missed 
the initial rapid acceleration phase of the CME. 

5. We estimate the rapid acceleration phase of the fil- 
ament between 17:50 UT and 22.33 UT, using the 
procedure described in section 3.2.f. This phase 
is marked by two dashed lines in the left panel of 
Figure [7] We fit functional forms of a parabolic, 
exponential and a cubic function to the true height- 
time profile of the filament apex during the rapid 
acceleration phase. These fits to the true height- 
time curves derived from two independent recon- 
struction methods described in section 3.2 and 3.3, 
respectively, are shown in the Figure [9] 

We now compare our analysis of the eruption kine- 
matics with the predictions of theoretical eruption mod- 
els described in the Introduction. We like to note that 
conclusions obtained from such a comparison should be 
read with some care and not be understood as a way to 
strictly confirm or rule out certain models. First, height- 
time data obtained with current instruments are still not 
accurate enough and typically do not have sufficient ca- 
dence to allow to pin down clearly the functional forms of 
rise profiles, which may behave very similar over the rel- 
atively short time scales of the initial rapid acceleration 
in solar eruptions. Also, a clear functional dependence 
may not be present if several acceleration mechanisms 
are at work simultaneously in an eruption. Second, for 
many models a proper investigation of the functional de- 
pendence of the eruption kinematics has not yet been 
reported, and even for most of those for which it was, 
there exists no parametric study, which may reveal kine- 
matics of a different functional dependence than reported 
for specific settings of the model parameters. 

Our data indicate that the rapid acceleration phase 
already started before the eruption became visible in 
the coronagraph data, so we restricted our analysis of 
the early acceleration phase to the filament observations. 
Our fits suggest that the filament enters an exponential 
rise phase at about 17:50 UT, which then appears to 
saturate from around 22:33 UT. Such exponential initial 
acceleration is line with many previous studies (see the 
Introduction) and supports the current picture that both 
quiescent and active region filament eruptions, and their 
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associated CMEs, are driven by the same mechanisms. 
It suggests the occurrence of an ideal MHD instability, 
here most likely the torus instability. We did not find 
indications of a clear writhing motion of the filament 
that would suggest the additional occurrence of the he- 
lical kink instability, although the twisted appearance of 
the filament sheet may indicate some untwisting of the 
magnetic field during the early phase of the eruption. 

The exponential accele ration found he r e is d ifferent 
from the cases studied by iSchriiver et al.1 ()2008j ) , where 
a cubic (or near-cubic) rise was found for two active re- 
gion filament eruptions. However, these authors showed, 
using numerical simulations, that a relatively large ini- 
tial velocity of the erupting structure at the onset of its 
rapid acceleration can change the subsequent rise behav- 
ior from exponential to cubic. The slow rise velocity of 
the filament (estimated from Figure [71 plateau in the fil- 
ament speed curve before the first vertical dashed line) 
before the filament enters rapid acceleration phase, is 
relatively small, about ~2. 5 km/s, however, com para- 
ble to the case described in ISchriiver et al.l (|2008h . The 
exponential rise al s o diff ers from the recent results by 
Uoshi fc Srivastaval (|2011l ) , who found a constant acceler- 
ation for both the slow rise and rapid acceleration phases 
of the two 3D-reconstructcd quiescent prominence erup- 
tions. However, these authors apparently did not fit func- 
tions other than parabolic, also the quality of their fits 
is not reported. 

While the data we considered here support the torus 
instability as the mechanism responsible for the initial 
rapid acceleration of the filament, they do not provide 
reasonable clues for the cause its preceding relatively 
long slow rise phase. We did not find indications of pre- 
flare brightening which are often used to draw conclu- 
sions. Hence, we do not find support for tether-cutting 
or magnetic breakout, but we cannot rule out the oc- 
currence of these and other reconnection-relatcd mecha- 
nisms, since the magnetic fields in the source region of 
the eruptions might have been simply too weak to pro- 
duce detectable brightening. We therefore refrain from 
speculating on the exact underlying mechanism respon- 
sible in the present case. However, more studies using 
the methods developed in this work and encompassing 
larger sets of observations, including the high-resolution 
SDO/AIA observations, could provide better clues. 
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TABLE 1 

Propagation Direction using Simple Triangulation and Marinus Method 
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a Anglcs a and /3 measured by simple triangulation method. 
b Carrington latitude and longitude using Marinus method. 
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Fig. 1.— The top panel shows the He 304 A filtergram observed from STEREO-B at 19:17 UT on 26 September 2009. LF and SF mark 
the large and the small filament (highlighted inside yellow rectangle, outlined by white line segments for clarity), respectively. The bottom 
panel shows a map of the radial magnetic field component, reconstructed from a synoptic MDI magnctogram, as STEREO-B would have 
seen it at the same time (Courtesy of Z. Mikic). The dashed lines indicate the polarity inversion line above which the filaments are located. 
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Fig. 2.— Four different views on the erupting filament in He II 304 A . The top panel shows the STEREO-Behind (left) and STEREO- 
Ahead (right) views. The bottom panel shows the Earth views by SoHO/EIT (left) and CORONAS/TESIS (right). 




Fig. 3. — STEREO-A (top) and B (bottom) images of the eruptive filament (LF and SF) are shown. The edges of the filament sheet 
(LF) arc outlined by dashed (red) and dotted (yellow) lines. The sheet appears to be twisted along its legs. 
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STEREO-B 




STEREO-A 



EARTH VIEW 
(EIT.TESISorSDO) 

Fig. 4. — An illustration of the de-projection method. The coronal structure is represented by a yellow loop extending above the solar 
limb. The top part of this loop is at an angle a (/}) to the plane-of-sky in the reference frame of STEREO-A (Earth View, EV). S is the 
separation angle between STEREO-A and EV. The segments ha and hb are the projected distances (R + /i)cos(/3) and (R + h)cos(a), 
respectively. S is related to a and by S = a + /J (if the loop is seen in front of the limb in one view and behind the limb in the other) or 
by S = | a — /3\ (if the loop is seen on the same side of the limb, i.e., either in front of or behind the limb, in both views). If ha, hb, and S 
are known, a and /3 can be determined. 



STEREO-A 19:16:15 UT 



SOHO EIT 19:19:38 UT 




Fig. 5. — Near-simultaneous STEREO-A and SoHO/EIT filtcrgrams in epipolar geometry. The segments h a and hf, are measured from 
the apex of the filament. The separation angle between Earth and STEREO-A is S = a + = 61°. 
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Fig. 6. — The top panels show the POS height-time profiles of the filament apex (left) and CME leading edge (right) as seen from 
different viewpoints (see inset for the observing satellites). The lower panels show the true height-time profiles derived by using the simple 
trigonometric method (Scction3.2). 
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Fig. 7. — The left (right) panels show, from top to bottom, the height, speed, acceleration and jerk (rate of change of acceleration) 
profiles of the filament apex (CME leading edge), respectively. In the top panels, the data points are shown by '+' marks and the solid 
line corresponds to a fourth order polynomial fit (H(t) = a + bt + ct 2 + dt s + et 4 ) to the data. The fitted profile H{t) is used to derive 
the speed, acceleration and jerk curves in the subsequent panels. The two vertical dotted lines in left panels correspond to the estimated 
duration of the rapid acceleration phase of the filament, which is fitted with different functional forms in the Figure [9] 
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Fig. 8. — The top, middle, and bottom panels show the three views of the sun from the STEREO-A, TESIS and STEREO-B, respectively 
in Marinus projection. The top part of the filament is marked by a white box in all images. 
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METHOD I - Quick Estimation 



METHOD II - 3D Reconstruction 
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Fig. 9. — The panels on the left and right show the altitude-time curve for the erupting filament derived from simple triangulation 
method (Section 3.2) and 3-D reconstruction method (Section 3.3). The four data points in red color in top right panel correspond to 
the altitude reconstructed by using SCC MEASURE method. The altitude-time curve correspond to the rapid acceleration phase of the 
filament eruption and is fitted for three functional forms viz. parabolic, exponential and cubic (from top to bottom). The reduced chi-square 
value of the fit is displayed on the top left corner of each panel. 



